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ABSTRACT : 


To investigate the effect of the role played by planetary beta in the generation and 
maintenance of the Leeuwin Current System (LCS), several numerical experiments are 
conducted with varying Coriolis parameterizations (plane or ®-plane). The seasonal 
effects of thermal gradients, wind forcing and North West Shelf (NWS) waters are also 
systematically explored in the numerical experiments. Additionally, the generation 
mechanisms for undercurrents along both the western and southern coasts of Australia are 
investigated. 

Although surface currents, undercurrents, and mesoscale activity are present in 
both the f-plane and §-plane experiments, those run on a $-plane show a more realistic 
depiction of these features due to the significant role played by the planetary beta effect in 
the LCS. The combination of thermal forcing, wind forcing, and NWS waters were found 
to be essential to maintain the highly seasonally varying LCS. Alongshore temperature 
gradients, which are generated at deeper levels by the strong flow of the Leeuwin Current, 
are sufficient to establish and maintain an equatorward (westward) undercurrent along 


western (southern) Australia. 
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I. INTRODUCTION 


The Leeuwin Current is a thermally driven, anomalous, surface eastern boundary 
current (EBC). It flows poleward over the continental shelf break off the western 
Australian coast, down to Cape Leeuwin, and then eastward into the Great Australian 
Bight (e.g., Cresswell and Golding, 1980). There is general agreement (e.g., Godfrey and 
Ridgway, 1985) that the Leeuwin Current is generated by a meridional pressure gradient 
resulting from the large amount of heating in the equatorial region and excessive cooling 
near the poleward region, which overwhelms the equatorward wind stress. The source 
for the Leeuwin Current is predominantly geostrophic inflow from the west (e.g., 
McCreary et al., 1986; Thompson, 1987) and is augmented by a souce from the North 
West Shelf (e.g., Gentilli, 1972), possibly having its origin in the Pacific Ocean (e.g., 
Hirst and Godfrey, 1993). As a result of this strong inflow of subtropical water towards 
the coast, the Leeuwin Current intensifies poleward (e.g., Batteen et al., 1992). The 
Leeuwin Current flows mainly but not exclusively in the austral autumn and austral 
winter (Church et al., 1989). The surface poleward current along the western Australian 
coast is weakest and shallowest between November and January, followed by the period 
of strongest and deepest flow between March and May. 

The Leeuwin Current is unlike any other subtropical EBC. The major subtropical 
EBCs such as the California, Canary, Peru, and Benguela currents are characterized by 
climatologically weak (<10 cm/s), broad (~1000 km wide), cold surface flow towards the 


equator in the direction of the prevailing winds, a poleward undercurrent, a shallow (<30 








m depth) thermocline, and high biological production due to vast regional upwelling 
(Parrish et al., 1983). Observational studies along the coast of Western Australia have 
shown that the Leeuwin Current is characterized by a strong (>150 cm/s at times), narrow 
(<100 km wide), poleward surface current that flows opposite the prevailing wind 
direction (e.g., Cresswell and Golding, 1980; Godfrey et al., 1986), anomalous warm 
water at the surface, a deep (>50 m depth) thermocline (Thompson, 1984), and lower 
biological production due to vast regions of downwelling (Batteen et al., 1992). 

Godfrey and Ridgway (1985) have examined the observed seasonal signals of the 
forcing mechanisms between 22°S and 32°S. They discuss for this region, the relative 
strengths of the equatorward wind forcing and the opposing pressure gradient (Figure 1), 
noting that the individual seasonal cycle of one forcing mechanism reinforces the 
seasonal trends of the other. Specifically, the pressure gradient term, which experiences a 
larger seasonal amplitude than wind forcing, reaches its strongest poleward forcing when 
the opposing wind stress reaches its weakest forcing. This combination of the two 
similar phases could explain the observed seasonal trends in the Leeuwin Current System 
(LCS). Added to this seasonal forcing is a surge of warm, less saline water which flows 
poleward from the North West Shelf (NWS) region during the austral autumn (1.e., March 
through May) (Smith et al., 1991). This surge of warmer water coincides with the period 
in which Godfrey and Ridgway (1985) found the strongest poleward pressure gradient 
forcing and weakest equatorward wind stress, resulting in a strong poleward driving 


function. 


Below the Leeuwin Current, there is an equally significant anomalous 
equatorward undercurrent off western Australia (e.g., Church et al., 1989). Smith et al., 
(1991) stated that the speed of the undercurrent can reach ~30 cm/s between ~250 m and 
350 m and observed that current meter data from the LUCIE experiment showed the 
equatorward current to be narrow and situated on the continental shelf slope between 
~250 m and 450 m. Although there is evidence of a westward flow of ~20 cm/s centered 
between ~400 m and 700 m depth at Cliffy Head off the southern coast of Australia 
(Cresswell and Peterson, 1993), no studies to date have clearly resolved whether there is 
a westward undercurrent in this region. | 

Previous numerical modeling studies by Batteen et al. (1992) investigated the 
effects of annual climatological wind forcing and initialized thermohaline gradients on 
the LCS, but the study was limited to the coast off western Australia, during the period 
when the LCS is observed to be strongest. Batteen and Butler (1998) examined the 
effects of continuously forced annual Indian Ocean thermohaline gradients on the LCS 
off both the western and southwestern coasts of Australia. 

This study extends prior efforts through the addition of seasonal (rather than 
annual) wind forcing and thermal gradients, permitting a more representative forcing of 
the highly seasonal LCS. In particular, a full primitive equation (PE) ocean model is 
used to investigate the individual and combined roles of seasonal wind stress and thermal 
gradients, including NWS forcing, in the generation of currents and eddies off the 
western and southern coasts of Australia. One particular focus of this study will be to 


examine the development and variation of the equatorward undercurrent along the 


western coast of Australia. The nature of the subsurface flow off the southwestern coast 
of Australia will also be investigated, with particular focus on whether there is a 
westward undercurrent. Additionally, the dependence of Coriolis parameterization (f- 
plane vs. f-plane) on the generation, evolution, and maintenance of currents and eddies 
off the western and southern coasts of Australia will be investigated. 

This study is organized as follows. In Chapter If we describe the numerical 


model and the specific experimental conditions. The results and energy analyses of the 


model experiments are presented in Chapter III. A summary is presented in section IV. 








II. MODEL DESCRIPTION 


A. MODEL EQUATIONS 


This study uses the eddy-resolving, limited-area EBC model of Batteen (1997) 
with open northern, western, and southern boundaries. Additionally, it has been modified 
by Batteen and Butler (1998) to have an open eastern boundary to account for interaction 
with the waters of the Great Australian Bight (GAB). 

The model is multi-level, uses non-adiabatic primitive equations on a B-plane or 
f-plane, and has both baroclinic and barotropic velocity components. It is based on the 
hydrostatic, Boussinesq, and rigid lid approximations. Equations governing the model 


are as follows: 
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In the above equations, ¢ is time, T is temperature, S is salinity, p is pressure, and p is 
density. A right-handed Cartesian coordinate system (x,y,z) is need where x points 
onshore, y alongshore, and z upward. The corresponding velocity components are (u,v,w). 
The variation of f with latitude is taken as f = (f, + By), where f, is the value of f at the 
mid-latitude of the region and 8B (or beta, the variation of the coriolis parameter with 
latitude) = Of / Gy is given the value at the mid-latitude of the model domain. Table 1 
provides a complete list of symbols found in the model equations, as well as values of 
constants used throughout the study. 

A space-staggered B-scheme (Arakawa and Lamb, 1977) is used for the 
horizontal finite differencing. Batteen and Han (1981) have shown this scheme to be 
appropriate when the grid spacing is approximately on the same order as, or less than, the 
Rossby radius of deformation. The horizontal grid spacing is 11 km in the east-west 
direction and 14 km in the north-south direction, while the internal Rossby radius of 
deformation is ~30 km, as determined by Feliks (1985). The model uses ten vertical 
layers, with constant z-levels, at depths of 10, 30, 75, 150, 250, 400, 600, 1226, 2283, and 
3656 m. This vertical spacing scheme concentrates more layers above the thermocline, 
which is in the dynamically active portion of the ocean. 

The model domain (Figure 2) encompasses the western and southwestern coasts 
of Australia, from 22.5°S to 40°S (1792 km alongshore), and from 107.5°E to 122.5°E 
(1408 km cross-shore). The coastal boundaries of the model domain are closed, and have 
both the tangential and normal components of velocity set to zero. Bottom topography 


has been omitted to focus on the roles of climatological wind forcing, thermal gradients, 





and Coriolis parameterization. The coastal boundary is modeled as a vertical wall and 
the constant depth used in the model is 4500 m. 

A modified version of the radiation boundary conditions of Camerlengo and 
O’Brien (1980) is used for the open boundaries to the north, south, east, and west with 
some spatial smoothing applied within five grid points (~60 km) of these boundaries 


(Batteen, 1997). Biharmonic lateral heat and momentum diffusion is used in the model 


17 om“/s) as in Batteen (1997). 


with the same choice of coefficients (i.e., 2.0 x10 
Holland (1978) showed that highly scale-selective biharmonic diffusion acts 
predominantly on submesoscales, while Holland and Batteen (1986) found that baroclinic 
mesoscale processes could be damped by Laplacian lateral heat diffusion. As a result, the 


use of biharmonic lateral diffusion should allow mesoscale eddy generation via 


barotropic (horizontal shear) and/or baroclinic (vertical shear) instability mechanisms. 


Weak (0.5 cm/s) vertical eddy viscosities and conductivities are also used. Bottom 
stress is parameterized by a simplified quadratic drag law (Weatherly, 1972), as in 
Batteen (1997). 

The method of solution is straightforward with the rigid lid and flat bottom 
assumptions because the vertically integrated horizontal velocity is subsequently 
nondivergent. The vertical mean flow can be described by a streamfunction which can be 
predicted from the vorticity equation, while the vertical shear currents can be predicted 
after the vertical mean flow is subtracted from the original equations. The other 
variables, i.e. temperature, salinity, density, vertical velocity, and pressure, can be 


explicitly obtained from the thermodynamic energy equation (6), salinity equation (7), 





equation of state (5), continuity equation (3), and hydrostatic equation (4), respectively 


(For more complete details on the method of solution, see Batteen, 1997). 
B. FORCING CONDITIONS AND EXPERIMENTAL DESIGN 


The model is forced from rest with climatological wind fields from a 2.5° by 2.5° 
grid of the European Centre for Medium Range Weather Forecasts (ECMWF) near- 
surface wind analyses (Trenberth et al., 1990). The monthly mean stresses based on 
twice daily wind analyses from 1980-1989 have been interpolated spatially to the 11 by 
14 km model resolution and temporally to daily wind values. Examples of the wind 
fields are shown in Figure 3, which illustrates the seasonal influence and migration of 
flow around the Southern Indian Ocean Subtropical High over the area encompassed by 
the model domain. In February (Figure 3a), the winds circulate anticyclonically around 
the Subtropical High positioned west-southwest of Cape Leeuwin. This flow produces 
upwelling favorable winds along both the western and southern coasts of Australia. By 
May (Figure 3b), the Subtropical High has migrated equatorward to ~30°S. The 
equatorward flow has weakened off western Australia and the wind has shifted to 
downwelling favorable westerlies off the southern coast of Australia. The Subtropical 
High reaches its maximum equatorward extent by July (Figure 3c), resulting in weak and 
variable flow off Western Australia. Strong westerlies and northwesterlies continue to 
dominate along the southern coast. By November (Figure 3d), the Subtropical High has 


migrated poleward, producing upwelling favorable winds with a strong equatorward 


component along the western coast of Australia. The winds off southern Australia have 
shifted from downwelling favorable westerlies to weak upwelling favorable easterlies. 

The effects of thermal gradients on the Leeuwin Current are incorporated by 
using the monthly temperature climatology from Levitus and Boyer (1994) to initialize 
the model, and every third day, to force the model at the western boundary (107.5°E). 
The seasonal temperature forcing conditions from Levitus and Boyer, which are initially 
assumed to be homogeneous, are shown for the seven levels in Figure 4a for the 
equatorward (22.5°S) and poleward (40°S) boundaries of the model domain. Note that at 
the poleward end of the western boundary, there are times (e.g., see months 4, 6, and 10 
in Figure 4a), when cooler water overlies warmer water in the upper layers. Before 
inserting such temperatures in the model, the temperatures have been slightly adjusted 
(e.g., Figure 4b) so that warm water always overlies the cool water. This is an 
initialization requirement in the model, since at a boundary there is only forcing, and no 
convective adjustment. As a result, the western boundary temperature forcing used in the 
model is that depicted in Figure 4b. 

Since the lower three levels of temperature from Levitus and Boyer (1994) do not 
exhibit much Borizenial variation, they are assumed to be constant for each level. The 
temperature values used for levels 8-10 are 3.88°C, 2.20°C, and 2.00°C. To focus on the 
role played by thermal forcing, salinity variations are omitted in this study. The constant 
salinity used is 34.7. 

While the temperatures to the north are warmer than those to the south, only the 


upper layer (e.g., Figure 4) temperature forcing conditions show significant seasonal 





variability with a temperature maximum in late summer to early fall and minimum in late 
winter to early spring througout the region. Below ~75m depth (not shown), both the 
seasonal temperature fluctuations and the temperature gradient weaken, as expected. 

The NWS waters provide a significant influx of warm tropical water to the 
Leeuwin Current and are used as an additional thermal forcing source. The climatology 
of Levitus and Boyer (1994) provides a representative temperature profile for the NWS 
waters, which is used in this study (Table 2). The NWS waters are modeled as a 
horizontal homogeneous raft of warm water simulating a “dam breaking” applied to the 
northeastern ocean boundary of the model domain once during each model year (see 
Figure 28b, for example). This model approach is consistent with the hypothesis of flow 
reversals as a function of monsoonal winds isolating water on the NWS, which 
subsequently flow poleward (Gentilli, 1972). 

Table 3 summarizes the design of the seven basic experiments and two sensitivity 
experiments. In Experiment 1, the effects of seasonal thermal gradients on the generation 
of the Leeuwin Current when the model is run on a f-plane 1s investigated. With the use 
of a numerical model, we examine in an idealized way, the development of flow from 
rest when the temperatures are initially assigned the offshore alongshore gradient. The 
system is then forced every third day at the western boundary (107.5°E) with seasonally- 
varying climatological temperature fields from Levitus and Boyer (1994). Experiment 2, 
a sensitivity experiment, is the same as Experiment | except an ideal (straightened) west 
and south coast are employed (e.g., see Figure 6a). To isolate the effects of planetary 


beta, Experiment 3 is run using the same initialization and forcing parameters as in 


10 





Experiment 1, but on a 8-plane. Experiment 4, a sensitivity experiment, is the same as 
Experiment 3 except only the top five layers are forced. Experiment 5 examines the role 
of seasonal wind forcing in the absence of thermal gradients on a f-plane. The model is 
initialized using a horizontally homogeneous temperature field and is then forced daily 
with seasonal ECMWF winds. Experiment 6 is the same as Experiment 5 except it is run 
on a B-plane. Experiment 7 investigates the combined effects of seasonal wind forcing 
and thermal gradients on a f-plane. Model initialization and forcing is a combination of 
Experiments 1 and 5. To isolate the effects of planetary beta, Experiment 8 is run using 
the same initialization and forcing parameters as in Experiment 7, but on a #-plane. 
Lastly, Experiment 9 examines the additional impact of the NWS waters (Table 2), along 
with seasonal wind forcing and thermal gradients. The model is run on a £-plane as in 
Experiment 8 with the release of NWS waters each year on model day 75 (which 
corresponds to mid-March, the timing of which is based on observations by Cresswell 
and Peterson, 1993). Note that in all of the experiments, the model is allowed to reach a 


quasi-equilibrium state by running each experiment for three years. 
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Ill. RESULTS FROM MODEL SIMULATIONS 


A. EXPERIMENTS 1 AND 2 —- RESULTS OF /-PLANE THERMAL FORCING 


EXPERIMENTS 


Due to the initial alongshore temperature field, the resulting pressure gradient 
establishes an onshore geostrophic inflow from the interior of the ocean. The onshore 
flow varies between ~2 and 10 cm/s. As the flow approaches the coast, it turns and forms 
a narrow poleward boundary current, which subsequently enters the Great Australian 
Bight (Figure 5a). Additional saeds flow augments the coastal current causing the 
magnitude of the current to increase in both the poleward and eastward directions (not 
shown). By day 27 (Figure 5b), maximum velocities of the order of ~65 cm/s are 
observed off Western Australia near Cape Leeuwin (see Figure 2 for geographical 
locations), and maximum velocities of ~100 cm/s are observed off southern Australia. 

The initial temperature field changes, as warm water is advected poleward and 
eastward throughout the domain. By day 27 (Figure 5b), 22°C water has moved as far 
south as Freemantle while 20°C water has been advected around Cape Leeuwin and 
reached the Great Australian Bight. Additionally, anticyclonic meanders develop along 
the coastal aeiiaiadia mesoscale features in the current intensify and grow, and by 
day 60 (see Figure 5c), eddies have formed off Freemantle, Cape Naturaliste and 


southwest of Cliffy Head. 
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Since this experiment was conducted on an f-plane, it was not clear that eddies 
would develop. In Batteen (1997), several experiments were run on an f-plane for the 
California Current System (CCS), the classical eastern boundary current system. In all 
experiments straight coastline geometry was used, while the model was forced using 
different types of wind stress. Under these conditions, eddies did not develop due to the 
undercurrent being too deep and weak, resulting in no near-surface baroclinic (vertical 
shear) or barotropic (horizontal shear) instability. Why did the eddies develop in this 
experiment? 

Since the eddies are generated at preferred locations (coastal promontories), is 
irregular coastline geometry the dominant mechanism for eddy formation? Experiment 2, 
which used a straightened coast (e.g., see Figure 6a), was conducted to ascertain whether 
eddy generation could occur when no irregular coastline geometry is used. As expected, 
eddies did form at the coastal promontory (see Figure 6a). In time (e.g., see day 159 in 
Figure 6b), there was additional eddy generation with no apparent geographical 
preference, which showed that eddies could occur regardless of the coastline geometry 
used. 

A meridional cross-section at ~32°S (see Figure 7) reveals a poleward surface 
current with a core velocity of ~50 cm/s overlying an equatorward undercurrent below 
~200 m depth with a core velocity of ~25 cm/s centered at ~500 m depth. The presence 
of a surface poleward current overlying a relatively shallow undercurrent results in both 


baroclinic and barotropic shear, which could lead to the subsequent generation of eddies. 
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What is the reason for the relatively shallow undercurrent in this case, since previous 
studies of the CCS did not show this feature? 

The formation of the equatorward undercurrent could possibly be attributed to a 
temperature gradient reversal with depth, i.e., Figure 4 shows that below ~400 m depth 
the temperature in the poleward end of the model domain is warmer than the temperature 
in the equatorward end of the model domain. This reversal in temperature gradient with 
depth offers a plausible explanation for the formation of the undercurrent (This topic will | 


be addressed in more detail in section N1.B.2). 


B. EXPERIMENTS 3 AND 4 - RESULTS OF 8-PLANE THERMAL FORCING 


EXPERIMENTS 
1. Spin-up Phase 


The 8-plane model simulations of the Leeuwin Current are quite similar to the f- 
plane runs, with anticyclonic meander and eddy formations at saustal promontories. As 
expected, the significant difference is the westward propagation of the eddies in the B- 
plane cases. For example, a comparison of Experiments 1 and 3 on day 87 (Figure 5d 
‘al Figure 8a, respectively) and on day 117 (Figure 5e and Figure 8b, respectively) 
shows the core of the eddy moving to the south at ~29°S in Experiment 1 (compare 
Figures 5d and 5e), while Experiment 3 shows the eddy west of Freemantle (originally at 


~111°E) propagating to the northwest at ~7 km/day (compare Figures 8a and 8b). 
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Additionally, westward propagation occurs along the southern coast of Australia in the 
8-plane experiment but not in the f-plane run (not shown). 

For the same reason that the Leeuwin Current intensifies as it flows poleward 
(geostrophic inflow), one would expect the undercurrent to intensify as it flows 
equatorward. Closer examination of cross-sections along the west coast however, reveals 
the opposite. For example, a comparison of day 15 for Experiment 3 at ~32°S (see Figure 
Qa) and at ~26°S (see Figure 9b) shows core velocities for the undercurrent decreasing 
equatorward from ~25 cm/s at ~32°S to 15 cm/s at ~26°S. A comparison for day 15 at 
the same latitudes for the f-plane Experiment 1 also shows the same trend (e.g., see 


Figures 10a and 10b). Why ts this? 


2. Analysis of Undercurrent Generation Mechanisms 


Previous modeling studies (Weaver and Middleton, 1988; Batteen eft al., 1992; 
Batteen and Butler, 1998) have generated an equatorward undercurrent off the western 
coast of Australia, consistent with observations of LUCIE (Boland et al., 1988). Off the 
southern coast of Australia, although Cresswell and Peterson (1993) have found evidence 
of a westward flow at Cliffy Head between 400 and 700 m depth with a speed of ~ 20 
cm/s, a lack of observations off this coast, particularly east of Cliffy Head, leaves some 
question as to the nature of the subsurface flow in that region, including its generation 


mechanism if present. 
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It is not clear what the mechanism for the generation of undercurrents is in the 
LCS. One possible mechanism for the development of the eguatorwatd undercurrent off 
western Australia could be the pressure gradient forcing at depth. As discussed earlier, 
the Indian Ocean seasonal thermal gradients experience a reversal between ~400 m and 
700 m depth (see Figure 4b), such that warmer waters are present in the poleward end of 
the model domain and cooler waters in the equatorward end. 

Another possible mechanism for the generation of the equatorward undercurrent 
is the following: When the system is in equilibrium, the strong flow of the surface 
current generates alongshore temperature gradients at deeper levels, with warm waters at 
the poleward end of the model domain off Cape Leeuwin and cooler waters at the 
equatorward end of the model domain. These gradients are sufficient to establish and 
maintain an equatorward undercurrent along Western Australia. 

An examination of the subsurface flow off the western coast of Australia in 
Experiment 3 shows that an equatorward undercurrent at ~ 400 to 600 m depth develops 
in the opposite direction of the sities current in this region. For example, Figure lla 
shows the equatorward flow at ~600 m depth. 

Of the two proposed mechanisms, the alongshore pressure gradient at depth 
cannot explain the equatorward undercurrent. In particular, the geostrophic velocity in 
response to the pressure gradient at depth would be in the offshore direction rather than 
onshore. Therefore we hypothesize that the second proposed mechanism, the generation 
of deeper alongshore thermal gradients by the surface flow, may explain the equatorward 


undercurrent. 
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To examine this further, we conducted a special experiment (Experiment 4 in 
Table 3) that removes the thermal gradient below 400 m depth. The model is initialized 
as in Experiment 3; however, the thermal forcing is only applied at the upper five levels. 
Level 6 (400 m) is given a constant temperature of ~10.5°C, while level 7 (600 m) is 
given the constant temperature of ~8.7°C. 

Figures 12a and 12b show cross-shore velocity sections for Experiment 4 at 
~32°S and at ~26°S for model day 3. At this time, the shallowest equatorward 
undercurrent with a core velocity of ~4 cm/s at ~600 m depth has developed near ~26° S. 
The undercurrent continues to intensify so that by day 9 (Figure 12c and 12d), 
equatorward undercurrents are well developed at both latitudes. The core of the deepest 
and strongest undercurrent (Figure 12c), however, is located near the poleward end of the 
model domain. The progression of development of this undercurrent in the absence of an 
initial pressure gradient, with the deepest and strongest current in the poleward region of 
the model domain, is consistent with our second proposed mechanism. In particular, the 
strengthening poleward surface flow generates alongshore temperature gradients at 
deeper levels. These temperature gradients, which are most pronounced in the region 
near Cape Leeuwin, are sufficient to establish and maintain an equatorward undercurrent 
off western Australia. 

We now examine the nature of the subsurface flow off the southern coast of 
Australia. Analysis of Figure 1la shows that a continuous westward undercurrent 


develops in the opposite direction of the surface flow off southern Australia in 


Experiment 3. A cross-shore velocity section at 117°E (Figure 13a) shows that the 
westward undercurrent has its core centered at ~600 m depth and has speeds of ~20 cm/s. 

Analysis of Figure 11b shows that there is also a continuous westward 
undercurrent off southern Australia in Experiment 4. A cross-shore velocity section at 
the same location (1.e., 117°E) for Experiment 4 (e.g., Figure 13b) also shows similar 
core depths and speeds as in Experiment 3. 

In Experiment 4, even with the absence of forcing of a deep thermal gradient, a 
westward undercurrent still develops. Can we attribute the development of this westward 
undercurrent along the southern coast of htiiiite to the generation of alongshore 
temperature gradients at deeper levels, as was seen in the development of the 
equatorward undercurrent along the western coast of Australia? 

Cross-shore velocity sections for Experiment 4 at 117°E (Figure 13b) and at 
122°E (Figure 13c) show core velocities for the undercurrent decreasing eastward from 
~15 cm/s at 117°E to ~7.5 cm/s at 122°E. Additionally, the core of the undercurrent 
deepens eastward from ~600 m depth at 117°E to ~1300 m depth at 122°E. Note that the 
deepest and strongest undercurrent in the westward region off southern Australia 
coincides with the region of the strongest surface eastward flow, i.e., Figure 13b (Figure 
13c) shows surface flows of ~60 cm/s (~40 cm/s). A comparison of velocity (Figures 
13b and 13c) and temperature (Figures 14a and 14b) cross-sections shows that the core of 
the surface current (undercurrent) at each location coincides with the downward (upward) 
bending of isotherms. These results are consistent with the proposed generation 


mechanism that alongshore temperature gradients, which are generated at deeper levels 
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by the strong flow of the Leeuwin Current, are sufficient to drive an alongshore flow and 


set up a westward undercurrent below ~400 m depth off southern Australia. 
C. EXPERIMENT 5 — RESULTS OF f-PLANE WIND FORCING EXPERIMENT 


In this f-plane experiment the model was run using seasonal wind forcing only 
(see Table 3), so that no thermal gradients were incorporated. While wind forcing only is 
an unlikely scenario for the LCS, it could still play an important role in the seasonal cycle 
of the LCS. 

As expected, the initially strong upwelling favorable wind stress results in an 
equatorward surface current off western Australia and a westward surface current off 
southern Australia (e.g., Figure 15a). A comparison of cross-sections of meridional 
velocity at ~32°S (Figure 16a) and ~26°S (Figure 16b) shows a shallow equatorward 
surface current above ~200 m depth with velocities decreasing equatorward from ~20 
cm/s at ~32°S, to ~10 cm/s at ~26°S. Beneath the surface current lies a much weaker 
poleward undercurrent with an average core depth of ~400 to 600 m and a speed of ~1 
cm/s, which is more developed near the equatorward end of the model domain (e.g., 
compare Figures 16a and 16b). Additionally, during the period of strong equatorward 
winds (for example, see Figure 3a), offshore Ekman transport has resulted in coastal 
upwelling off the western coast (e.g., Figure 15a). The weak easterlies off the southern 


coast have also produced colder water consistent with upwelling along with a westward 


20 








current (e.g., Figure 15a), which attains a maximum velocity of ~50 cm/s near Cliffy 
Head (not shown). 

As the westward surface current off the southern coast flows around Cape 
Leeuwin, turns equatorward, and flows along the western Australia coast, it encounters 
the coastline northwest of Dongara, which deflects the surface flow to the northwest. 
This deflection of the flow sets up an anticyclonic meander (e.g., Figure 15b), and 
subsequently develops into a warm core, anticyclonic eddy at ~29°S (e.g., = Figure 
15c). Concurrently, a cold core, cyclonic eddy forms to the south of the warm core eddy, 
resulting in the formation of an eddy pair (e.g., see Figure 15c). | 

Another interesting development, similar to a study conducted by Batteen and 
Murray (2000) in which anticyclonic mesoscale activity occurred in an embayment, i.e., 
the Gulf of Cadiz, occurs in this case. By day 297 (see Figure 15d) an anticyclonic eddy 
has formed inshore of the equatorward surface current in the embayment area to the north 
of Cape Naturaliste and persists throughout the remainder of the year. The strong 
equatorward surface current just offshore and the lack of the planetary 8-effect (which 
would result in westward propagation of the eddy) essentially traps the eddy, preventing 
any offshore movement. With the exception of the eddy formations just discussed, there 
was no other mesoscale activity during the first year of the model run. 

It is interesting to note that only a relatively weak subsurface current below ~1000 
m depth develops off the southern coast. For example, Figure 17a shows that the 
subsurface flow on day 243 is found below the westward surface current of ~30 cm/s at 


~1200 m depth, and is relatively weak, with a typical eastward speed of ~4 cm/s. By the 
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end of the year (e.g., Figure 17b), the equatorward flow below the deep surface current of 
~65 cm/s has attained a maximum speed of ~8 cm/s, which is relatively weak compared 
to the intense surface current. 

By the end of the year (e.g., Figure 15e), due to the return of strong equatorward 
winds, upwelling occurs everywhere inshore of the coastal current off western Australia 
with the most intense upwelling occurring equatorward of ~31°S. Additionally, the 


equatorward surface current has intensified to speeds as strong as ~80 cm/s. 


D. EXPERIMENT 6 —- RESULTS OF 8-PLANE WIND FORCING EXPERIMENT 


As discussed earlier, results from previous sensitivity studies along eastern ocean 
boundaries (e.g., Batteen, 1997) have produced various results depending on the type of 
wind forcing, coastline geometry and Coriolis parameterization used. For the CCS, 
Batteen (1997) forced the model from rest using seasonal climatological winds along a 
straight California coast on both an f-plane and $-plane. Her study also included 
seasonal climatological wind forcing on an irregular coastline using B-plane only and 
focused on the role of the beta effect in generating, evolving and maintaining classical 
mesoscale features in the CCS. The results of her study showed that a surface 
equatorward current, upwelling of cooler water along the coast, and a poleward 
undercurrent develop. In addition, both baroclinic and barotropic instability develop due 
to the vertical and horizontal shear between the equatorward surface current and 


poleward flow, resulting in the generation of meanders near the coast, upwelling 


Ze 


filaments extending offshore, and the subsequent development of seasonal or semi- 
permanent cyclonic and anticyclonic eddies farther offshore. While the meridional 
variability of f (B-plane) and variability in the equatorward component of wind stress — 
were shown to be necessary ingredients for generating realistic vertical and horizontal 
structures for the cores of the surface equatorward and subsurface poleward currents, 
eee in the coastline geometry were shown to be important for “anchoring” 
upwelling and filaments, and for enhancing the growth of meanders and eddies. The 
results from several numerical experiments in this process-oriented study supported the 
hypothesis that wind forcing and coastline irregularities on a $-plane are important 
mechanisms for generation of observed mesoscale features of the CCS (Batteen, 1997). 

The expected outcome of wind forcing only on a B-plane in our experiment in the 
LCS would be similar to what Batteen (1997) found in the CCS. Additionally, one would 
not expect to see the development of as strong of an equatorward surface current (which 
formed in the fplane run) along the western coast of Australia, due to the westward 
dispersion of energy associated with the planetary beta effect. 

As expected, in response to the strong upwelling favorable wind stress, an 
equatorward surface current off western Australia and a westward surface current off 
southern Australia (e.g., Figure 18a) develop, similar to the f-plane case. A comparison 
of meridional cross-sections for day 33 at ~32°S and ~26°S (see Figures 19a and 19b, 
respectively) reveals weaker equatorward surface current and stronger poleward 
undercurrent speeds than in the previous experiment (e.g., see Figures 16a and 16b). This 


is due to the beta effect, which results in weaker surface currents and stronger 
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undercurrents. Other differences from the f-plane are seen with the development of 
cyclonic meanders off Cape Leeuwin, off the promontory nonhiweet of Freemantle, and 
off the promontory northwest of Dongara (e.g., Figure 18b). 

By day 165 (see Figure 18c), during the period of weakest winds off western 
Australia (e.g., see Figure 3c), the equatorward flow has weakened significantly from ~15 
cm/s to ~7.5 cm/s (e.g., compare Figures 20a and 20b). While the westward surface 
current off the southern coast has also weakened, the onshore Ekman transport produced 
by the stronger westerlies in that region (e.g., see Figure 3c) continues to augment the 
westward flow, extending the current west of Cape Leeuwin (e.g., Figure 18c). Asa 
result, a large (~250 km diameter) cyclonic eddy has formed off Cape Leeuwin and the 
surface flow “jets” westward and equatorward around this feature. Additionally, 
anticyclonic eddies develop resulting in the formation of eddy pairs (e.g., see eddy pairs 
at ~28°S and at ~30°S in Figure 18c). 

With the return of upwelling favorable equatorward winds (e.g., see Figure 3d) 
during the latter part of the year along the western Australia coast, a strong equatorward 
surface current and a poleward undercurrent re-develop. For example, cross-shore 
velocity sections at ~26°S (e.g., Figure 20c) show a relatively strong (~30 cm/s) 
equatorward surface current overlying a weaker (~15 cm/s) poleward undercurrent 
centered at ~400m depth. 

The structure of the westward surface current and eastward undercurrent along the 


southern coast of Australia are similar to that of the /-plane experiment, 1.e., a moderate 
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surface current continues to overly a weaker undercurrent. This is expected since the 


planetary beta effect should have no effect on a southern coast, which acts like an f-plane. 


E. EXPERIMENT 7 - RESULTS OF f-PLANE THERMAL AND WIND 


FORCING EXPERIMENT 


In previous experiments, the presence of seasonal thermal gradients only caused 
geostrophic inflow to dominate throughout the domain. The presence of wind forcing 
only, due to the seasonally varying alongshore and cross-shore components of wind 
velocity (Figure 3), produces a more complex flow regime. In particular, at the 
equatorward end of the model domain, since the winds have offshore and equatorward 
components, Ekman theory predicts poleward and offshore transport of mass, 
respectively. At the poleward end of the model domain off Western Australia, since the 
winds have onshore and equatorward components, there should be equatorward and 
offshore transport, respectively. Off the southern coast, the change from upwelling 
favorable to downwelling favorable winds as the eae change should produce offshore 
and then onshore mass transport. 

The combination of seasonal Indian Ocean temperature forcing at the western 
boundary with seasonal wind forcing is expected to produce results similar to available 
observations. Based on the results of the previous experiments, it is speculated that the 
seasonally-varying strengths of the thermal gradients and wind forcing mechanisms (see 


Figure 1) will result in an unstable poleward flow that adjusts to seasonal variations. 
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As expected, the initial thermal gradient, with temperature decreasing polewards, 
establishes a predominantly onshore geostrophic inflow in the interior ocean. On 
approaching the coast, the onshore flow turns and forms a poleward boundary current 
(e.g., Figure 21a), which advects warmer water from the equatorward end of the model 
domain. Offshore to the west, the equatorward wind stress dominates and drives offshore 
flow. Along the southwest coast, the poleward surface current rounds Cape Leeuwin and 
flows eastward towards the Great Australian Bight. The upwelling favorable winds east 
of Albany (see Figure 3a), initially drive an offshore flow, while the westerlies poleward 
of 38°S drive an onshore flow. 

By day 30 (Figure 21b), the poleward flow is well established and intensifies 
poleward. For example, a cross-shore section of meridional velocity near Freemantle 
(Figure 22a), depicts the relatively narrow (~20 km) structure of the poleward current 
extending to ~250 m depth with a core surface velocity of ~40 cm/s. There is also a 
weaker (~25 cm/s) equatorward undercurrent centered at ~SOO m depth. A comparison 
of this cross-section with one at ~26°S (Figure 22b) reveals a weaker poleward surface 
current (~4 cm/s) to the north, as expected when thermal forcing dominates. 

The current “jets” around Cliffy Head, accelerating toward the Great Australian 
Bight region (e.g., Figure 21b). Warm (20°C) water has been advected around Cape 
Leeuwin and reached Albany. A cross-shore section off the southern coast (Figure 23) 
depicts an eastward surface current, which has a core velocity of ~60 cm/s, overlying a 


westward undercurrent, which has a core velocity of ~15 cm/s. 
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By day 30 (Figure 21b), anticyclonic meanders appear, near Freemantle, and ‘in 
the vicinity of Cape Naturaliste, Cape Leeuwin, and Cliffy Head. In time, the 
anticyclonic meanders develop into eddies (e.g., see eddy at ~31°S in Figure 21c) and 
propagate slightly offshore (e.g., see eddy at ~31°S in Figure 21d), to a similar offshore 
extent as was seen in Experiment 1, which had thermal forcing only. 

By the end of the year, wind forcing has become the dominant forcing mechanism 
over the entire mode! domain. A cold, westward “jet” has formed and rounded Cape 
Leeuwin (e.g., Figure 21e), similar to Experiment 5 (see Table 3), and extends up the 
coast to Freemantle. A cross-shore section of meridional velocity at ~26°S (see Figure 
24a) reveals that the flows have reversed (e.g., compare with Figure 22b), resulting in an 
intense (~55 cm/s) and deep (extends to ~1100 m depth) equatorward surface current 
overlying a weaker (~6 cm/s) poleward undercurrent. Once established, the intense 
equatorward “jet” that forms in the f-plane experiment due to the wind forcing dominates 
over any poleward surface flow that could be due to the effects of thermal forcing (e.g., 


years 2 and 3, not shown). 


F. EXPERIMENT 8 ~ RESULTS OF §$-PLANE THERMAL AND WIND 


FORCING EXPERIMENT 


In time, the addition of the planetary beta effect should result in a more realistic 
depiction of eddies and coastal upwelling off the western coast of Australia due to the 


westward dispersion of energy. Initially, however, the ®-plane model simulation is quite 
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similar to the fplane run of Experiment 7 (see Table 3), with anticyclonic meander and 
eddy formations (e.g., compare Figures 21b and 25a of Experiments 7 and 8, 
respectively) at similar geographical locations (for example, off coastal promontories). A 
comparison of meridional cross-sections at ~32°S and ~26°S (e.g., compare Figures 26a 
and 26b, respectively) also reveals comparable poleward surface current and equatorward 
undercurrent speeds and core depths as seen at both locations in the previous experiment 
(e.g., see Figures 22a and 22b). 

Similar to the thermal forcing results of Experiments 3 and 4, a significant 
difference from the f-plane run is the westward propagation of the eddies in the ®-plane 
case. For example, a comparison of Experiments 7 and 8 on day 63 (Figure 21c and 
Figure 25b, respectively) and on day 123 (Figure 21d and Figure 25c) shows the core of 
an eddy at ~31°S meandering slightly to the north in Experiment 7 (e.g., compare Figures 
21c and 21d), while Experiment 8 shows a similarly located eddy propagating to the 
west-northwest at ~5 km/day (e.g., compare Figures 25b and 25c). 

By day 165 (Figure 25d), during the period when the Leeuwin Current is observed 
to be the strongest, a poleward flow is maintained at ~25-45 cm/s off the western coast. 
A cross-shore velocity section near Freemantle (Figure 24b) depicts the typical nature of 
the Leeuwin Current. A shallow (< 250 m depth) poleward surface current with a speed 
of ~25 cm/s, which is confined to within 50 km of the coast, overlies a weaker 
equatorward undercurrent with a speed of ~15 cm/s at ~500 m depth. 

By the end of the model year (e.g., Figure 25e), when the winds have become 


strongly equatorward off the western coast, the signature of the surface poleward 
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Leeuwin Current is lost, and equatorward currents dominate the coastal region. A 
comparison of meridional cross-sections at ~32°S and ~26°S (see Figures 27a and 27b, 
respectively) depicts the dominance of the wind forcing. At ~32°S, a shallow (< 200 m 
depth) equatorward surface flow of ~15 cm/s, sia a poleward undercurrent with a 
core speed of ~10 cm/s at ~400 m depth is evident, while at ~26°S a poleward 
undercurrent with a core speed of ~7.5 cm/s has shoaled inshore of a broad equatorward 
surface flow of ~15 cm/s (Figure 27b). This is due to the planetary beta effect, which has 
acted to move water offshore along the west ad of Australia. In conjunction with the 
equatorward surface current, coastal upwelling along the northwest coast (e.g., north of 
~27°S) is more prevalent in this case than in the f-plane experiment (e.g., compare 
Figures 21e and 25e). 

Further examination of years 2 and 3 (not shown) reveal a very weak Leeuwin 
Current signature for only a short period of time during the austral winter, a period when 
it should be most pronounced. Otherwise, wind forcing dominates, resulting in an 
unrealistic depiction of the LCS similar to the results of the previous experiment. Based 
on these results, we speculate that there must be another forcing mechanism that, in 
conjunction with thermal forcing, acts on the LCS and overcomes the dominant 


equatorward currents due to the wind forcing. 
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G. EXPERIMENT 9 —- RESULTS OF 8-PLANE THERMAL, WIND, AND 


NORTHWEST SHELF WATERS FORCING EXPERIMENT 
1. Spin-up Phase 


The introduction of warm NWS waters, in conjunction with wind forcing and 
thermal gradient forcing, is expected to generate the most realistic depiction of the 
Leeuwin Current System. The release of the NWS waters during mid-March of each 
model year is expected to introduce baroclinicity into the Leeuwin Current (e.g., Batteen 
et al., 1992), although the effects should diminish away from the source region. Batteen 
and Butler (1998) demonstrated that the addition of warmer NWS waters could intensify 
the current velocity and temperature fronts, and act to produce a more unstable flow 
regime. 

Before the onset of NWS waters, the flow is, as expected, the same as in 
Experiment 8 (e.g., see Figure 25a), with onshore geostrophic flow driving a poleward 
surface current along the western coast (not shown). Off southern Australia, the current 
“jets” around Cape Leeuwin and flows into the Great Australian Bight. By day 75 (not 
shown), anticyclonic meanders and instabilities have developed at coastal promontories, 
anticyclonic eddies have propagated westward advecting warm water to the offshore 
regions, and the flow around the southern coast has weakened in response to easterly 


winds. 
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The release of the NWS waters in the coastal region north of Carnarvon at day 78 
(Figure 28a) introduces a strong thermal front in the region and au gments the poleward 
flow. By day 102 (Figure 28b), the “raft” of NWS water has been advected poleward and 
has reached Dongara. A cross-shore section of velocity at ~26°S (e.g., Figure 29) depicts 
the well-developed surface current with speeds of ~60 cm/s overlying a weak (~6 cm/s) 
undercurrent centered at ~250 m depth. The flow near Cape Leeuwin has become 
stronger than in Experiment 8, attaining a speed of ~60 cm/s (e.g., see Figure 28b), 
consistent with the increased thermal forcing and the weakening equatorward winds 
during this time. 

The effects of the NWS waters have reached their farthest extent by day 165 
(Figure 28c). A comparison of Experiments 8 and 9 depicts a more well-developed 
Leeuwin Current along the western and southern coasts of Australia in Experiment 9 than 
in Experiment 8 (e.g., compare Figures 28c and 25d, respectively). In Experiment 9 
(Figure 28c), the warm surge of NWS water is contained close to the coast with 24°C 
water extending south of Dongara and a pocket of 22°C water approaching Cliffy Head, 
while Experiment 8 (Figure 25d) results show 20°C water at both of these locations. 
Cross-sections of meridional velocity (Figure 30) at 30°S show a poleward surface flow 
of ~60 cm/s for Experiment 9 (Figure 30a) while the cross-section for Experiment 8 
(Figure 30b) depicts a weaker (~25 cm/s) poleward flow. Cross-sections of zonal 
velocity (Figure 31) at 117°E show an easterly surface flow of ~120 cm/s and a westerly 


undercurrent of ~30-40 cm/s for Experiment 9 (Figure 31a), while the cross-section for 
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Experiment 8 (Figure 31b) shows a weaker (~70 cm/s) easterly surface flow overlying a 
stronger westerly undercurrent of ~30-40 cm/s. 

By day 363 (Figure 28d), when the opposing wind forcing (see Figure 3) has 
returned, the effects of the NWS waters diminish and the poleward surface current off the 
western coast is no longer discernible. A comparison of meridional cross-sections at 
~32°S (Figure 32a) and ~26°S (Figure 32b) depicts the equatorward surface flow along 
the western coast of Australia with a stronger (~30 cm/s) equatorward surface current 
near Shark Bay (Figure 32a) than at Freemantle (~10 cm/s) (Figure 32b) due to stronger 
equatorward wind forcing (e.g., see Figure 3d). The opposing wind off the southern coast 
also results in the weakening and subsequent reversal of the surface flow from eastward 


to westward (e.g., Figure 28d). 
2. Quasi-equilibrium Phase 


A continuation of the model run shows that a quasi-steady state has been reached 
by year 3, and reveals that most features recur on a seasonal basis. Specifically, the 
release of NWS waters on day 807 (day 78 of year 3 , not shown) tends to augment the 
thermal] forcing, which results in the development of a relatively strong surface poleward 
current. Anticyclonic meanders form near Shark Bay and Dongara and subsequently 
develop into eddies (e.g., Figure 33a). By day 915 (day 185 of year 3, Figure 33b), 
during the period when the Leeuwin Current is observed to be the strongest, the water off 


the coast of Dongara is 22°C, consistent with observations by Cresswell and Peterson 
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(1993), while 19°C water has been advected around Cape Leeuwin, extending to the east 
of Cliffy Head. In addition, a typical meridional cross-section during the same period at 
~30°S (e.g., see Figure 34) depicts a 25 cm/s poleward surface current inshore of a 15 
cm/s core equatorward undercurrent centered at ~250 m depth along the western coast, 
comparable to results from the LUCIE experiment which measured 38.5 cm/s for the 
poleward surface current and 5.9 cm/s for the undercurrent with a core located at ~300 m 
depth during the austral winter (e.g., Church et al., 1989; Boland et al., 1988; Smith et al., 
1991). A typical zonal cross-section along the southern coast at 117°E (e.g., see Figure 
35a) reveals an ~80 cm/s eastward surface current, consistent with observations reported 
by Cresswell and Peterson (1993), overlying a deeper, 10 cm/s westward undercurrent. 

In time, as expected, the effects of the NWS water begin to diminish and by late 
spring-early summer, the return of opposing winds (to the flow of the Leeuwin Current) 
diminishes the poleward and eastward flow. Subsequently, the flow reverses and eddy 
generation decreases (e.g., Figure 33c), pre-existing eddies propagate offshore (see eddy 
at 34°S in Figures 33c and 33d), and upwelling occurs along the entire coast (e.g, Figures 
33c and 33d), and is strongest north of 28°S. This scenario completes the seasonal cycle 
of the Leeuwin Current System, as the wind forcing causes the Leeuwin current to 
weaken and retreat equatorward. The strengthening of the thermal gradient the following 
year, in conjunction with the release of NWS waters in mid-March, causes the Leeuwin 
Current to re-intensify and advance poleward during the austral autumn and summer. 

Cross-sections of meridional velocity (Figure 36) at 30°S show a weak poleward 


undercurrent of ~4 cm/s for Experiment 9 (e.g., Figure 36a) while the cross-section for 
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Experiment 6 (e.g., Figure 36b) depicts a stronger (~15 cm/s) poleward undercurrent. 
During this period when wind forcing is, once again, dominant along the west coast of 
Australia, the influence of the opposing surface forcing mechanisms due to thermal 
forcing and NWS waters is still apparent, resulting in a weaker poleward undercurrent on 
the west coast in Experiment 9. 

Cross-sections of zonal velocity at 117°E show a relatively strong eastward 
undercurrent of ~40 cm/s for Experiment 9 (e.g., Figure 35c), while the cross-section 
(e. g., Figure 37) for Experiment 6 depicts a much weaker (~10 cm/s) and deeper eastward 
undercurrent. Along the southern Australia coast, the westward surface flow in both 
Experiments 6 and 9 has developed in response to strong westward winds, accompanied 


by an opposing eastward undercurrent which develops below the westward flow. 


3. Analyses of Mean and Eddy Kinetic Energy 


Typical cross-sections for the surface poleward and subsurface equatorward 
currents along the west coast (e.g., Figure 34) show considerable horizontal and vertical 
shear in the upper layer currents. To determine the type of instability that could generate 
the meander and eddy features, more detailed analysis was performed. Barotropic 
instability can result from horizontal shear in the currents, while vertical shear ig the 
currents can cause baroclinic instability. 

To indicate locations of barotropic and baroclinic instability, horizontal maps of 


the upper layer mean kinetic energy (MKE), eddy kinetic energy (EKE), and temperature 
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contours with velocity vectors were time-averaged for the six month period when the 
Leeuwin Current was the strongest (April-September). MKE and EKE results for this 
period (e.g., Figures 38a and 38b) show that high values of MKE and EKE are found 
along the jet axis of the surface current off the western and southern coasts of Australia, 
particularly in the regions of eddy development (e.g., compare Figures 38a and 38b to 
Figure 39a). High values of EKE are also indicated offshore, where eddies propagate. A 
comparison of Figures 38a and 38b shows that maximum values of MKE and EKE occur 
in the same general region. This is consistent with the results of Batteen (1997), which 
showed that eddies in the CCS are generated from instabilities of the mean surface 
current and the undercurrent via barotropic and/or baroclinic instability processes. 
Additionally, the NWS waters have added greater horizontal and vertical shear, resulting 
in greater barotropic and baroclinic instability processes, all along the west coast. 

To help elucidate the forcing mechanism(s) responsible for the MKE and EKE 
areas (Figure 38) depicted in icaiaiaian 9, horizontal maps of the upper layer MKE, 
EKE, and temperature contours with velocity vectors, time averaged for the six month 
period from April — September, were also obtained for Experiments 3, 6, and 8. In the 
thermal forcing only case (Experiment 3), typical sreaessecaoi of meridional velocity 
(e.g., Figure 40a) at ~30°S and zonal velocity (e.g., Figure 40b) at ~117°E both show 
horizontal and vertical shear in the upper layers near the coast. The MKE and EKE 
results (e.g., Figure 41) show that the highest MKE and EKE values are found along the 


jet axis of the eastward surface current (e.g., see Figure 39b) along the southwest and 
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southern coasts of Australia, particularly in the regions near Cape Leeuwin and Cliffy 
Head. | 

In the wind forcing only case (Experiment 6), a typical cross-section of 
meridional velocity (e.g., Figure 42a) at ~30°S shows horizontal shear in the upper layers, 
while a typical cross-section (e.g., Figure 42b) of zonal velocity at ~117°E shows no 
near-surface horizontal or vertical shear. The MKE results (e.g., Figure 43a) show that 
the highest values are found along the jet axis of the westward surface current off the 
south coast of Australia, with another, less intense maximum near Dongara (e.g., see 
Figure 39c). The highest values of EKE (e.g., Figure 43b) are located off the western 
coast of Australia between ~29°S to ~27°S and offshore of Capes Leeuwin and 
Naturaliste (e.g., Figure 39c). The absence of EKE (1.e., see Figure 43b), in conjunction 
with no near-surface shear (e.g., Figure 42b), along the southern coast of Australia in 
Experiment 6 indicates that the eddy development along the southern coast of Australia 
in Experiment 9 must be due to thermal rather than to wind forcing. 

In Experiment 8 (thermal and wind forcing), typical cross-sections of meridional 
velocity (e.g., Figure 44a) at ~30°S and zonal velocity (e.g., Figure 44b) at ~117°E show 
both horizontal and vertical shear in the upper layers. The MKE results (e.g., Figure 45a) 
show that the highest values are found along the jet axis of the eastward surface current 
along the south coast of Australia and along the jet axis of the equatorward surface 
current offshore of Shark Bay and between ~33°S to ~25°S west of ~112°E (e.g., compare 
Figure 45a to Figure 39d). The EKE results (e.g., Figure 45b) show highest values in the 


same regions as the MKE results, but depict an eddy development region all along 
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western Australia. Comparison of EKE results from Experiment 8 (e. g., Figure 45b) and 
Experiment 9 (e.g., Figure 38b) reveal higher values of EKE all along the west coast in 
Experiment 9, due to the addition of NWS waters. The NWS waters have added greater 
horizontal and vertical shear, resulting in greater barotropic and baroclinic instability 


processes, all along the west coast. 
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IV. SUMMARY 


The objective of this study was to use a multilevel PE model to investigate the 
effect of seasonal thermal gradients and wind forcing and the role played by planetary 
beta in the generation and maintenance of the LCS. Toward this end, several numerical 
experiments were run with varying Coriolis parameterizations (f-plane or B-plane) using 
either thermal gradients only, wind forcing only, or a combination of both. Additionally, 
the seasonal effects of the NWS waters were included in another experiment, yielding 
insight into the importance of this secondary forcing mechanism. A final goal of this 
research was to investigate the development of undercurrents along the western and 
southern coast of Australia. 

The results of Experiment 1 revealed how seasonal thermal gradients on an f- 
plane generate and maintain a poleward and eastward surface current overlying a 
westward and equatorward undercurrent. Additionally, eddies were generated at 
preferred locations (coastal promontories) but did not migrate westward. To determine 
whether irregular coastline geometry was the dominant mechanism for eddy formation, 
Experiment 2 was run using a straightened coast. Eddy generation still occurred, but with 
no apparent geographical preference. The formation of eddies was attributed to the 
baroclinic and barotropic aie between the thermally forced poleward surface current 
and equatorward undercurrent. 

In Experiment 3, the beta effect, as expected, did not yield any major differences 


from Experiment 1, other than the westward propagation of eddies. Examination of 
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meridional cross-sections along the west coast revealed that the undercurrent weakened 
equatorward (opposite of the poleward intensifying surface current due to geostrophic 
inflow) when one would have expected the undercurrent to intensify due to geostrophic 
inflow. This result led to an investigation of possible generation mechanisms for the 
coastal undercurrent. A comparison of model results from Experiment 4, which, after 
inititialization used no thermal forcing below the upper five levels, with Experiment 3 
showed that along the western and southern coasts of Australia, the strong flow of the 
surface current generated alongshore temperature gradients at deeper levels, which were 
sufficient to establish and maintain an equatorward (westward) undercurrent along 
western ection Australia. 

The results from Experiment 6, the beta plane, wind forcing only case, 
demonstrated the importance of spatially and temporally varying wind in the seasonal 
nature of the LCS. The strong equatorward winds of the austral spring and summer 
generated significant upwelling along the entire coastline, particularly in the equatorward 
region of the model domain. This upwelling became weak and somewhat suppressed 
with the weaker winds of the austral fall and winter, which corresponds to the period 
when the Leeuwin Current is observed to be the strongest. Cyclonic meander formation 
and subsequent eddy formation off Cape Leeuwin and the major promontories were 
consistent with previous EBC studies (e.g., Batteen, 1997). 

To isolate the role of the beta effect, the results of Experiment 5 (/-plane) were 
compared to Experiment 6. The fplane case also generated a westward flow along the 


southern coast and equatorward flow along the west coast. Closer examination, however, 
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revealed a stronger equatorward surface current and a weaker poleward undercurrent than 
the ®-plane case. This was attributed to the beta effect, which results in weaker surface 
currents and stronger undercurrents due to offshore transport at the surface and the 
resultant inflow at depth. As the model year progressed, an intense equatorward “jet” 
developed, which was stronger and farther offshore than in Experiment 6. An interesting 
feature of the f-plane case was the development and continued presence of an 
anticyclonic eddy inshore of the “jet” in the embayment north of Cape Naturaliste. With 
the exception of the eddy just discussed and an eddy pair that generated off the 
promontory northwest of Dongara, mesoscale activity was non-existant in the f-plane 
case. 

The results of Experiment 7 explored the combined role of seasonal thermal 
gradients and wind forcing, on an f-plane, in the LCS. The relative strength of the 
combined forcing mechanism initially varied with the seasonal fluctuations of thermal 
gradients and wind forcing. The poleward and eastward Leeuwin Current strengthened 
and weakened in response to these seasonal variations. Anticyclonic meanders were 
generated and, in time, developed into eddies which propagated to a similar offshore 
extent as was seen in Experiment 1, the thermal forcing only, f plane case. By the end of 
the year, when wind forcing became dominant over the entire model domain, a cold, 
westward “jet” developed along the south coast and moved around the cape, extending up 
the coast to Freemantle. Once this “jet” formed, it did not allow for the establishment of 
any sustained poleward surface flow, even when thermal forcing was the dominant 


forcing mechanism. 
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The addition of the beta effect in Experiment 8 provided improved results over 
Experiment 7. The early spin-up phase was quite similar to the Polane case, with the 
exception of the westward propagation of eddies. A cold, coastal “jet” did not form and 
the addition of the beta effect allowed for more realistic upwelling to occur at the end of 
the year, when surface flow became equatorward along the west coast due to the 
dominance of the wind forcing. In the austral fall and winter of years 2 and 3, however, 
the Leeuwin Current was only weakly depicted and of short duration when it should have 
been most pronounced. 

The added annual contribution of NWS waters in Experiment 9 augmented the 
onshore geostrophic inflow, strengthening the poleward flow and advecting warmer 
waters along the coastal region. The warmer waters along the coast resulted in advection 
of warmer waters offshore by anticyclonic eddies. The steadily increasing forcing from 
stronger thermal gradients and the poleward surge of NWS waters combined to extend 
the duration of the strong poleward and eastward surface flow of the Leeuwin Current 
through the austral winter, which recurred annually, resulting in the most realistic 
depiction of the LCS. Analysis of the surface current and undercurrent revealed realistic 
speeds and a similar core depth for the undercurrent when compared to actual 
observations. Additionally, the development of the undercurrent along the west and 
south Australia coasts was consistent with the generation mechanisms investigated in 
Experiment 4. Specifically, alongshore temperature gradients, which are generated at 
deeper levels by the strong flow of the Leeuwin Current, are sufficient to establish and 


maintain an equatorward (westward) undercurrent along western (southern) Australia. 
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Analyses of MKE and EKE was performed to determine the type of instability 
that could generate the meander and eddy features in the various experiments. In 
Experiment 9, a comparison of MKE and EKE revealed that the highest values of both 
occurred in the same general region, mainly along the jet axis of the surface current off 
the western and southern coast of Australia. The eddies were found to be generated from 
instabilities of the mean surface current and the undercurrent via barotropic and/or 
baroclinic instability processes, consistent with the cess of Batteen (1997). 

This study provided insight into the interaction of seasonal thermal gradients and 
wind forcing in the LCS off the western and southwestern coast of Australia, as well as 
the role of planetary beta in the generation and maintenance of the LCS. The model 
qualitatively reproduced the complex circulation of this region, while yielding insight 


into the role of each forcing mechanism as well as identifying the important seasonal 


effects in the LCS. 
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Figure 1. Annual cycle of net forcing for the Leeuwin Current between 22°S and 32°S: 
Dashed line is observed pressure gradient forcing; dotted line is observed wind stress; 
solid line is net forcing. Positive is southward (from Godfrey and Ridgway). 
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Figure 2. The model domain for the Leeuwin Current System (LCS) is bounded by 


22.5°S to 40°S, 107.5°E to 122.5°E. The cross-shore (alongshore) resolution is 11 km (14 
km). Geographical locations and prominent features are labeled. 
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Figure 3. Climatological (1980-1989) ECMWF winds in m/s for: (a) February, (b) May, 
(c) July, and (d) November. Maximum wind vector is 10 m/s. 
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Figure 4. Plots of (a) actual and (b) slightly adjusted Levitus climatological temperature 


fields used as seasonal forcing along the western boundary (107.5°E) at 22.5°S (top) and 
40°S (bottom). 


51 





Temperature Forcing Along Western Boundary (107.5 E) 





N 
o1 


N 
© 


Temperature (C) 


Temperature (C) 





Months (January — December) 


(b) 


52 








DEPTH (m) : 10 
T (DAY) : 9 





26.0°S 
30.0°S 
Lid 
= 
— 
= 
;— 
a 
54.0°S 
38.0°S 
uO 9 ba 112.0°E 114.0°EF 116.0°E 118.0°E 120.0°E 
iv 10, LONGITUDE 


Temperature (deg C) and Velocity (cm/s) 
(a) 


Figure 5. Temperature contours and velocity vectors at 10 m depth at days (a) 9, (b) 27, 

(c) 60, (d) 87, and (e) 117 of experiment 1. The contour interval is 1°C. To avoid clutter, 
the velocity vectors are plotted every third grid point in the east-west direction and every 
fourth grid point in the north-south direction. Maximum velocity is 100 cm/s. 
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Figure 6. Temperature contours and velocity vectors at 10 m depth at days (a) 51 and (b) 
159 of Experiment 2. The contour interval is 1°C. To avoid clutter, velocity vectors are 
plotted every third grid point in the east-west direction and every fourth grid point in the 
north-south direction. Maximum velocity is 100 cm/s. 
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Figure 7. Cross-shore section of meridional velocity (v) at ~32°S for Experiment 2 on 


day 33. Contour interval is 10 cm/s for poleward flow (dashed lines) and 5 cm/s for 
equatorward flow (solid lines). 
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Figure 8. Temperature contours and velocity vectors at 10 m depth at days (a) 87 and (b) 


117 of Experiment 3. The contour interval is 1°C. To avoid clutter, the velocity vectors 


are plotted every third grid point in the east-west direction and every fourth grid point in 
the north-south direction. Maximum velocity is 100 cm/s. 
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Figure 9. Cross-shore section of meridional velocity (v) at (a) ~32°S and (b) ~26°S for 


Experiment 3 on day 15. Contour interval is 10 cm/s for poleward flow (dashed lines) and 
5 cm/s for equatorward flow (solid lines). 
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Figure 10. Cross-shore section of meridional velocity (v) at (a) ~32°S and (b) ~26°S for 


Experiment 1 on day 15. Contour interval is 10 cm/s for poleward flow (dashed lines) and 
5 cm/s for equatorward flow (solid lines). 


65 





DEPTH (m) 





LATITUDE : 26.18 
T (DAY) & 15 


200 


400 


600 





800 


1000 


112.80°9E 113.00°E 113.20°E 113.40°E 113.60°E 113.80°F 
LONGITUDE 


Meridional Velocity (cm/s) 
(b) 


66 





114.00°E 








DEPTH (m) : 600 
T (DAY): 15 


26 OS ee OS at ee 


‘ 
a 
4 
4 
4 
4 
4 
a 
a 
a 








¥ ¥ * - » ~ 
30.073 + 
* a * 
—o . » * * » » » * e 
mar “ «< « a < « « “ «a a - * 
* - * rd ~ > * ~ * * » * 
+ v ’ * Aa 
34.0°S 
. - - -~ « 7 = = Sd bd Lal hd = 
a a * o rs « « « « * a « 
* od » > * » * ° ~ ~ cad * 
“a » » a » * . . « @ 4 
a a a a 4 a a a a a a a A a a a a 4 4 4 a a a a 4 4 4 4 a s 
* . a) * » » . * » > » » » » » * » » » . > » . ¥ ’ ? v » » * . 1d 
-« “a a = «a — = = = = = = = -= Ld a = = = = = =a = «a = a a a a = at =a a 
38 O°S ¥ v vr * * » » . . . + * + * o * o o a Po ro o o + * * * ra rd » Cd 


110: 0°E 192:0°E 114.0°E 116.0°E 113:0°E 120.0°E 
UVa = 50D LONGITUDE 


Velocity (cm/s) 
(a) 
Figure 11. Velocity vectors at 600 m depth for (a) Experiment 3 and (b) Experiment 4 at 
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Figure 12. Cross-shore section of meridional velocity (v) at (a and c) ~32°S and (b and d) 
~26°S for Experiment 4 on day (a and b) 3 and (c and d) 9. Contour interval is 5 cm/s for 
poleward flow (dashed lines) and (a and c) 2.5 cm/s and (b and d) 1 cm/s for equatorward 
flow (solid lines). 
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Figure 13. Cross-shore section of zonal velocity (u) at (a and b) 117°E and (c) 122°E for 


(a) Experiment 3 and (b and c) Experiment 4 on day 15. The contour interval is (a and b) 


5 cm/s and (c) 2.5 cm/s for westward flow (dashed lines) and 10 cm/s for eastward flow 


(solid lines). 
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Figure 14. Cross-shore section of temperature at (a) 117°E and (b) 122°E for Experiment 
4 on day 15. The contour interval is 1°C. 
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Figure 15. Temperature contours and velocity vectors at 10 m depth at days (a) 33, (b) 
195, (c) 231, (d) 297, and (e) 363 of Experiment 5. The contour interval is 1°C. To avoid 
clutter, the velocity vectors are plotted every third grid point in the east-west direction 
and every fourth grid point in the north-south direction. Maximum velocity is 100 cm/s. 
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Figure 16. Cross-shore section of meridional velocity (v) at (a) ~32°S and (b) ~26°S for 


Experiment 5 on day 33. Contour interval is 1 cm/s for poleward flow (dashed lines) and 
5 cm/s for equatorward flow (solid lines). 
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Figure 17. Cross-shore section of zonal velocity (u) at 117°E for Experiment 5 on days 
(a) 243 and (b) 363. Contour interval 1s 5 cm/s for westward flow (dashed lines) and 2 


cm/s for eastward flow (solid lines). 
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Figure 18. Temperature contours and velocity vectors at 10 m depth at days (a) 33, (b) 


90, and (c) 165 of Experiment 6. The contour interval is 1°C. Maximum velocity for (a) 
and (b) is 100 cm/s and for (c) is 50 cm/s. 
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Figure 19. Cross-shore section of meridional velocity (v) at (a) ~32°S and (b) ~26°S for 


Experiment 6 on day 33. Contour interval is 1 cm/s for poleward flow (dashed lines) and 
5 cm/s for equatorward flow (solid lines). 
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Figure 20. Cross-shore section of meridional velocity (v) at 34°S at days (a) 33 and (b) 
165 and at ~26°S at day (c) 315 for Experiment 6. Contour interval is (a) 1 cm/s, (b) 2.5 
cm/s, and (c) 5 cm/s for poleward flow (dashed lines) and (a and c) 5 cm/s and (b) 2.5 
cm/s for equatorward flow (solid lines). 
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Figure 21. Temperature contours and velocity vectors at 10 m depth at days (a) 15, (b) 


30, (c) 63, (d) 123, and (e) 363 of Experiment 7. The contour interval is 1°C for (a) — (d) 
and 2°C for (e). Maximum velocity is 100 cm/s. 
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Figure 22. Cross-shore section of meridional velocity (v) at (a) ~32°S and (b) ~26°S for 
Experiment 7 on day 30. Contour interval is (a) 10 cm/s and (b) 1 cm/s for poleward flow 
(dashed lines) and 5 cm/s for equatorward flow (solid lines). 
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Figure 23. Cross-shore section of zonal velocity (u) at 117°E for Experiment 7 on day 
30. Contour interval is 5 cm/s for westward flow (dashed lines) and 10 cm/s for eastward 
flow (solid lines). 
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Figure 24. Cross-shore section of meridional velocity (v) at (a) ~26°S for Experiment 7 
on day 345 and (b) ~32°S for Experiment 8 on day 165. Contour interval is (a) 1 cm/s and 
(b) 5 cm/s for poleward flow (dashed lines) and 5 cm/s for equatorward flow (solid lines). 
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Figure 25. Temperature contours and velocity vectors at 10 m depth at days (a) 30, (b) 

63, (c) 123, (d) 165, and (e) 363 of Experiment 8. The contour interval is 1°C for (a) — (d) © 
and 2°C for (e). Maximum velocity is 100 cm/s. Note: In (e), temperature contours 
intersecting coast north of 28°S are 16°C. 
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Figure 26. Cross-shore section of meridional velocity (v) at (a) ~32°S and (b) ~26°S for 


Experiment 8 on day 30. Contour interval is (a) 5 cm/s and (b) 2.5 cm/s for poleward 
flow (dashed lines) and (a) 5 cm/s and (b) 2.5 cm/s for equatorward flow (solid lines). 


110 











DEPTH (m) 


LATITUDE. : 26.18 


T (DAY) : 30 
= ~ ~ 
“ Xx ——— 
Tg ed ~ i x mM Teas 
~ bare ~—_ ~ ~e 
~ ~ bar ~ 
~e ae ~ ~ 
~ ~ ~ 
ee bears a = 
Sy ~ S ™ Ne 
~~ O.. Pa f 
~e Hh é 
~ ON ae é ‘ 
100 ae = ; F 











300 
200 
700 
900 
1100 ° 
tas a © ad 113.30°E 1135.50°E LiTL37 0" 113.90°E 
LONGITUDE 


Meridional Velocity (cm/s) 
(b) 


111 








Sac 


363 


LATITUDE : 
T (DAY) 


O 
O 
8) 


1000 








1200 


115.60°E 1715.80°E 


115.40°E 


LONGITUDE 
| Velocity (cm/s) 


T1S.00°E 115.20°E 


114.80°E 


lOna 


* 


Meri 


(a) 


Figure 27. Cross-shore section of meridional velocity (v) at (a) ~32°S and (b) ~26°S for 
Experiment 8 on day 363. Contour interval is 2.5 cm/s for poleward flow (dashed lines) 


and 5 cm/s for equatorward flow (solid lines). 


112 





26313 
363 


LATITUDE : 
T (DAY) 


70 


200 


“~~ 
~ 
~ 
“a 
~ 


- 
~— 
- 
se 
= a 
~ a 
ma 


¢7~ 
o” 
y 
“ 
\ 
1} 
\ \. 
\ 
~ 
oO © O O 
oO oO Oo oO 
t+ tO oO O 


- 

~~ 
———— 
——— 
ad 
-_——o 
wnwe7orn"™ 
= 
n= 
ad 
Pr aad 
o_o 
—— 





1200 


112:8°E 113 .2°E 11S66°E 114.0°E 
LONGITUDE 


112:4°E 


1A2Z:Or8 


(b) 
113 


Meridional Velocity (cm/s) 








DEPTH (m) : 10 
T (DAY): 78 


26.0°S 


30:0°S 


LATITUDE 


34.0°S 








110.0°E 112.0°E 114.0°E 116.0°E 118.0°E 120.0°E 
UV “100; LONGITUDE 


Temperature (deg C) and Velocity (cm/s) 
(a) 
Figure 28. Temperature contours and velocity vectors at 10 m depth at days (a) 78, (b) 


102, (c) 165, and (d) 363 of Experiment 9. The contour interval is 1°C for (a), (b), and (d) 
and 2°C for (c). Maximum velocity is 100 cm/s. 
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Figure 29. Cross-shore section of meridional velocity (v) at ~26°S for Experiment 9 on 
day 102. Contour interval is 10 cm/s for poleward flow (dashed lines) and 1 cm/s for 
equatorward flow (solid lines). 
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Figure 30. Cross-shore section of meridional velocity (v) at 30°S for (a) Experiment 9 
and (b) Experiment 8 at day 165. Contour interval is (a) 10 cm/s and (b) 5 cm/s for 
poleward flow (dashed lines) and 5 cm/s for equatorward flow (solid lines). 
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Figure 31. Cross-shore section of zonal velocity (u) at 117°E for (a) Experiment 9 and 
(b) Experiment 8 at day 165. Contour interval is 5 cm/s for westward flow (dashed lines) 
and 10 cm/s for eastward flow (solid lines). 
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Figure 32. Cross-shore section of meridional velocity (v) at (a) ~32°S and (b) ~26°S for 


Experiment 9 on day 363. Contour interval is 2.5 cm/s for poleward flow (dashed lines) 
and (a) 2.5 cm/s and (b) 5 cm/s for equatorward flow (solid lines). 
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Figure 33. Temperature contours and velocity vectors at 10 m depth at days (a) 870, (b) 


915, (c) 1026, and (d) 1092 of Experiment 9. The contour interval is 2°C for (a), 3°C for 
(b), and 1°C for (c) and (d). Maximum velocity is 100 cm/s. 
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Figure 34. Cross-shore section of meridional velocity (v) at 30°S on day 948 for 
Experiment 9. Contour interval is 5 cm/s for poleward flow (dashed lines) and 5 cm/s for 
equatorward flow (solid lines). 
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Figure 35. Cross-shore section of zonal velocity (u) at ~117°E for Experiment 9 on day 
(a) 930, (b) 990, and (c) 1074. Contour interval for velocities greater than (less than) 10 
cm/s is 10 cm/s (5 cm/s) for westward flow (dashed lines). For eastward flow(solid lines), 
the contour interval is 10 cm/s in (a) and (c) and 5 cm/s in (b). 
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Figure 36. Cross-shore section of meridional velocity (v) at 30°S for (a) Experiment 9 
and (b) Experiment 6 at day 1074. Contour interval is (a) 2 cm/s and (b) 5 cm/s for 
poleward flow (dashed lines) and 10 cm/s for equatorward flow (solid lines). 
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Figure 37. Cross-shore section of zonal velocity (u) at 117°E on day 1074 for 
Experiment 6. Contour interval is 10 cm/s for westward flow (dashed lines) and 2.5 cm/s 


for eastward flow (solid lines). 
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Figure 38. Horizontal map at 10 m depth of average (a) mean kinetic energy (MKE) and 
(b) eddie kinetic energy (EKE) for Experiment 9 in model year 3, time averaged for April 
through September, when the Leeuwin Current is strongest. Contour interval for kinetic 
energy ((cm/s)*) greater than 1000 is 500 and less than 1000 is 250. 
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Figure 39. Temperature contours and velocity vectors at 10 m depth for Experiment (a) 
9, (b) 3, (c) 6, and (d) 8 in model year 3, time averaged for April through September, 
when the Leeuwin Current is strongest. The contour interval is 1°C. Maximum velocity is 


100 cm/s. 


138 





DEPTH (m) : 10 





26.0°S 
30.0°S 
Lid 
(am 
— 
= 
— 
= 
34.0°S 
38.0°S 
oli 8 bag ot 112.0°E 114.0°E 116.0°E 118.0°E 120.0°E 
U_AVG_PERIOD , V_AVG_PERIOD ————>> 100. LONGITUDE 


Average Temperature (deg C) and Velocity (cm/s) 
(b) 


139 





DEPTH (m): 10 





26.0°S 

30.0°S 
Lit 
Cc 
=) 
= 
< 

34.0°S 

110.0°E 172.0°E 114.0°E 116.0°E 118.0°E 120.0°E 
U_AVG_PERIOD , V_AVG_PERIOD —————>> 100. LONGITUDE 


Average Temperature (deg C) and Velocity (cm/s) 
(c) 


140 





DEPTH (m) : 10 


LATITUDE 





110.0°E T1205 114.0°E 116.0°E 118.0°E 120.0°E 


U_AVG_PERIOD , V_AVG_PERIOD ————> 100. LONGITUDE 
Average Temperature (deg C) and Velocity (cm/s) 
(d) 


14] 





: 30S 
; 930 


LATITUDE 
T (DAY) 





200 


4 


{0 


400 


O 
O 
+6) 


1000 








1200 


114.4°E 114.8°E 115.2°E 


114.0°E 


LONGITUDE 


| Velocity (cm/s) 


(a) 


lOna 


* 


Meri 


Figure 40. Cross-shore section of (a) meridional velocity (v) at 30°S on day 930 and (b) 
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zonal velocity (u) at 117°E on day 918 for Experiment 3. Contour interval is 5 cm/s for 
(a) poleward and (b) westward flow (dashed lines) and 5 cm/s for (a)equatorward and (b) 


eastward flow (solid lines). 
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Figure 41. Horizontal map at 10 m depth of average (a) mean kinetic energy (MKE) and 
(b) eddie kinetic energy (EKE) for Experiment 3 in model year 3, time averaged for April 
through September, when the Leeuwin Current is strongest. Contour interval for kinetic 
energy ((cm/s)”) greater than 1000 is 500 and less than 1000 is 250. 
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Figure 42. Cross-shore section of (a) meridional velocity (v) at 30°S on day 912 and (b) 
zonal velocity (u) at 117°E on day 924 for Experiment 6. Contour interval is (a) 5 cm/s 


and (b) 10 cm/s for (a) poleward and (b) westward flow (dashed lines) and (a) 10 cm/s 
and (b) 2 cm/s for (a) equatorward and (b) eastward flow (solid lines). 
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Figure 43. Horizontal map at 10 m depth of average (a) mean kinetic energy (MKE) and 
(b) eddie kinetic energy (EKE) for Experiment 6 in model year 3, time averaged for April 
through September, when the Leeuwin Current is strongest. Contour interval for kinetic 
energy ((cm/s)*) greater than 1000 is 500 and less than 1000 is (a) 250 and (b) 125. 
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Figure 44. Cross-shore section of (a) meridional velocity (v) at 30°S on day 915 and (b) 
zonal velocity (wu) at 117°E on day 930 for Experiment 8. Contour interval is (a) 2.5 cm/s 
and (b) 10 cm/s for (a) poleward and (b) westward flow (dashed lines) and (a) 2.5 cm/s 
and (b) 10 cm/s for (a) equatorward and (b) eastward flow (solid lines). 
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Figure 45. Horizontal map at 10 m depth of average (a) mean kinetic energy (MKE) and 
(b) eddie kinetic energy (EKE) for Experiment 8 in model year 3, time averaged for April 
through September, when the Leeuwin Current is strongest. Contour interval for kinetic 
energy ((cm/s)”) greater than 1000 is 500 and greater than 1000 is 250. 
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Table 1. Values of constants used in the model 


[Constant] Value Definition SsS—SSS 
3 ___[347______| Constant Reference Salty 


7.5 x 107 Saline Expansion Coefficient 
Number of Levels In Vertical 
1.1x 10°cm Cross-Shore Grid Spacing 
1.4x 10°cm Alongshore Grid Spacing 
4.5 x 10°cm Total Ocean Depth 
-0.76 xX 107s" Mean Coriolis Parameter 
Acceleration of Gravit 
Biharmonic Momentum Diffusion Coefficient 
Biharmonic Heat Diffusion Coefficient 


Vertical Eddy Viscosity 
Vertical Eddy Conductivity 
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Table 2. Temperature profile for North West Shelf Waters 


Layer (Depth in meters) Temperature (°C) 
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Table 3. Summary of specific experimental design — 


7 Thermohaline | Wind _ | Ideal (Straight) | North West 
Parameterization Gradient Forcing Coastline Shelf 
Water 
Els 


1 

8-plane 

i 
(top 5 layers) 

| YES N 

N 

N 


O O 

S O 

O O 
ia : ° 
| 5 |  fplane YES 
i ae O NO 
eae O 

O NO 

O 



















[| B-plane | = =NO_ | YES [| =NO_ | NO 
7 | fplane | YES | YES | =NO | NO 
| YES | YES |  =NO | NO | 
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